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mixing  ratio  profiles.  This  improvement  enables  the  Kuo  scheme  to  parameterize  shallow 
to  medi-m  clouds  as  well  as  deep  clouds.  Tests  of  the  Kuo  scheme  using  a  semi-prognos¬ 
tic  approach  and  a  Cloud  Cluster xModel  indicate  that  the  improved  version  verifies 
better  with  observation  during  weaker  convective  periods. 

A  set  of ^experimental  predictionswere  made  and  the  results  clearly  demonstrated 
the  ability  of  the  AFGL  model  to  predict  large-scale  stratiform  precipitation.  With 
thej  ir.cor porat ion  of  the  ^present?’ modified  Kuo  scheme,  the  area  of  convective  precipi¬ 
tation  can  also  be  well  predicted.  However,  the  predicted  convective  precipitation 
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area  is  generally  broader  and  the  amount  of  rainfall  smaller  than f obser vat ionsy-  -Th-rs 
may  be  attributed  to  the  resolution  of  the  AFGL  model  in  that  it  cannot  resolve 
the  narrow  band  of  the  cold  front,  the  major  mesoscale  rain-producing  system.  The 

implementation  of  the  A-5  scheme  in  the  AFGL  model  also  produced  a  reasonable  distribu¬ 
tion  of  convective  precipitation^  The  precipitation  area  is  more  concentrated  and 
I  s-oeomes  it'  is  proaucec  m  an  observed  clear  area.  This  result  may^-have  beeh  caused 
|  by  t re  current  implementation  of  the  A-3  scheme  in  the  AFGL  model,  in  which  the  cloud 
base  is  assumed  to  be  5GC  m  above  the  ground  surface.  This  assumption  is  valid  ir.  the 
trcoics.  but  it  may  produce  unrealistic  convective  rainfall  in  cry  areas. 
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One  of  the  central  problems  in  numerical  weather  prediction  is  to 
relate  the  sub-grid  scale  condensation,  evaporation,  and  vertical  transports 
of  heat  and  moisture  by  cumulus  clouds  to  the  variables  predicted  by  the 
model.  This  problem  is  known  as  the  cumulus  parameterization.  The  most 
important  element  of  any  cumulus  parameterization  scheme  is  the  prediction 
of  the  rainfall  rate  and  the  vertical  distributions  of  the  convective  heat¬ 
ing  and  moistening  profiles.  The  rainfall  rate  is  directly  related  to  the 
total  amount  of  heat  added  to  and  moisture  removed  from  the  large-scale 
system  by  condensation  and  evaporation.  The  vertical  distributions  of  heat 
and  moisture  by  clouds  can  change  the  large-scale  stability  and  moisture 
structure  of  the  atmosphere  which,  in  turn,  can  affect  the  future  develop¬ 
ment  of  large-scale  weather  systems. 

Several  cumulus  parameter i zat ion  schemes  have  been  proposed  in  the 
pa3t.  The  simplest  and  most  commonly  used  scheme  i3  the  one  proposed  by  Kuo 
0965,  1974)  based  on  the  assumption  that  the  convective  precipitation  is 
related  to  the  vertically  integrated  moisture  convergence.  In  this  scheme, 
usually  only  one  type  of  cloud  is  produced  at  a  given  time  and  the  cloud 
heating  and  moistening  profiles  are  proportional  to  the  temperature  and 
humidity  differences  between  the  cloud  and  the  environment.  Another  scheme 
gaining  Increasing  attention  is  the  one  proposed  by  Arakawa  and  Schubert 
(1974)  which  takes  into  account  the  coexistence  of  a  spectrum  of  clouds. 
The  quasi-equilibrium  approximation,  which  requires  that  clouds  stabilize 
the  atmosphere  as  the  large-scale  motion  generates  moist  convective  insta¬ 
bility,  permits  the  determination  of  cloud  properties. 
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Parallel  to  the  effort  of  the  Air  Force  Geophysics  Laboratory  (AFGL) 
in  developing  a  Moist  Global  Forecast  Model,  we  carried  out  a  research  pro¬ 
ject  in  cumulus  parameterization  under  the  sponsorship  of  the  Meteorology 
Division  of  the  AFGL.  The  objectives  of  this  research  are:  (1)  to  criti¬ 
cally  evaluate  the  Kuo  and  Arakawa-Schubert  (A-S)  cumulus  parameterization 
schemes  for  numerical  weather  prediction;  and  (2)  to  improve  these  cumulus 
parameterization  schemes  for  the  purpose  of  improving  precipitation  fore¬ 
casts  on  the  global  scale.  The  results  of  our  research  on  the  A-S  scheme  of 
cumulus  parameter i zat ion ,  including  an  improved  computational  algorithm  for 
computing  the  cloud  spectrum  under  the  quasi-equilibrium  approximation,  have 
been  presented  in  Scientific  Report  No.  1  (Kao  and  Ogura,  1985).  Part  I  of 
this  report  contains  the  results  of  our  research  on  the  Kuo  scheme.  The 
major  improvements  to  the  Kuo  scheme  include  the  prediction  of  cloud  top 
height  and  the  incor porat ion  of  the  effect  of  entrainment  on  cloud  tempera¬ 
ture  and  mixing  ratio  profiles.  This  improvement  enables  the  Kuo  scheme  to 
parameterize  shallow  to  medium  clouds  as  well  as  deep  clouds.  Tests  of  the 
Kuo  scheme  using  a  serai-prognostic  approach  and  a  Cloud  Cluster  Model  indi¬ 
cate  that  the  improved  version  verifies  better  with  observation  during 
weaker  convective  periods.  Part  II  contains  the  results  of  applying  the 
improved  Kuo  scheme  and  the  A-S  scheme  of  cumulus  parameterization  in  the 
AFGL  Moist  Global  Forecast  Model.  In  contrast  to  the  forecast  using  the 
original  code  of  cumulus  parameter  i  zation  in  the  AFGL  model,  both  the  im¬ 
proved  Kuo  scheme  and  the  A-S  scheme  predicted  the  area  of  convective  pre¬ 
cipitation  reasonably  well. 
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PART  ONE 

A  Cumulus  Parameterization  Study 
with  Special  Attention  to  the  Kuo  Scheme 

1 •  Introduction 

It  is  recognized  that  the  latent  heat  released  in  cumulus  convection 
and  the  vertical  transports  of  heat  and  moisture  in  clouds  play  key  roles  in 
determining  the  structure  of  the  temperature  and  moisture  fields  in  the 
atmosphere.  However,  since  the  horizontal  scale  of  cumulus  clouds  is  at 
least  an  order  of  magnitude  smaller  than  the  grid  scale  used  in  representing 
large  and  mesoscale  flows,  the  net  effect  of  convection  and  condensation 
upon  the  large-scale  flows  must  be  parameter  i  zed  in  terms  of  the  known 
large-scale  variables. 

Several  cumulus  parameterization  schemes  are  used  in  large-scale  and 
mesoscale  models.  Each  scheme  is  different  in  a  way  directed  by  convenience 
and  by  the  type  of  problem  and  integration  model  involved.  Tne  most  impor¬ 
tant  elements  of  any  cumulus  parameterization  scheme  are  the  prediction  of 
the  rainfall  rate  and  the  convective  heating  and  moistening  profiles.  The 
rainfall  rate  determines  the  total  amount  of  heat  added  to  the  large-scale 
system  by  condensation,  while  the  vertical  distribution  of  heat  and  moisture 
by  clouds  can  change  the  large-scale  stability  and  moisture  structure  in 
ways  that  can  either  enhance  or  suppress  the  development  of  large-scale 
systems.  The  simplest  and  most  popular  scheme  of  cumulus  parameterization 
is  based  on  the  assumption  tha-  convective  preoipi tat  ion  is  related  to  the 
vertically  integrated  moisture  convergence.  Kuo  (1965,  197^)  parameter i zed 

the  release  of  latent  heat  by  relating  convective  activity  to  the  moisture 
converging  in^o  a  conditionally  unstable  atmospheric  column.  This  report 
will  focu3  on  this  scheme.  Section  2  reviews  some  of  tne  Kuo-type  eumul  is 


_  o  _ 


parameter i nation  schemes ,  while  Section  3  presents  a  semi-prognostic  test  of 
the  Kuo  scheme.  Section  M  describes  the  present  modified  Kuo  scheme,  an  : 
Section  5  presents  the  results  of  the  semi-prognostic  studies  with  the  pres¬ 
ent  scheme.  Section  6  presents  the  prognostic  model  results  using  a  flout 
'luster  Model. 

2.  Review  of  Cumulus  Parameterization  Schemes  Based  on  Moisture  Conver¬ 
gence 

A  cumulus  parameterization  based  on  the  total  column  moisture  conver¬ 
gence  was  first  proposed  by  Charney  and  Eliassen  (1964).  The  application  of 
this  method  to  tropical  cyclones  has  established  the  concept  of  CISK  (Condi¬ 
tional  Instability  of  the  Second  Kind).  In  this  method,  cumulus  clouds  are 
permitted  to  penetrate  deeply  into  the  unsaturated  atmosphere  in  areas  of 
low-level  mass  convergence.  Extensive  use  of  this  method  has  led  to  several 
related  but  not  identical  techniques.  Kuo  (1965,  1 97  ^  refined  this  concept 
by  distributing  cumulus  heat  and  moisture  vertically  with  a  simple  cloud 
mr-  VI .  The  amount  cf  cumulus  heating  is  assumed  to  be  proportional  t  o  the 
vertically  Integrated  moisture  convergence.  This  technique  has  often  been 
r-.  f  erred  to  as  trie  Kuo  scheme  of  cumulus  parameterization.  In  this  section, 
we  will  review  several  variations  of  Kuo's  parameterization. 

2.1.  Kuo  (1965): 

Kuo's  suo-cme  stipul  it.en  tha*  cumulus  convection  always  o : curs  in  re- 
gi.ns  of  deep  layers  of  conditionally  unst  able  stratification  and  mean  low- 
lev-  oonvergeri,'e  ,r  <  ■  M  i‘  the  vertic.il  profiles  of  temperature  and  w.»‘. 

vapor  mixing  ratio  inside  the  cloud  follow  those  cf  a  non-entraining  moist 
ill  it  i'll’  process  .  Trie  .••_-rti  c  il  i  1  s’ r  i  but  i  or...  of  con  ve  ;t  1  ve  heat  ing  uni 
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moistening  can  be  determined  by  the  mixing  of  environmental  and  cloudy  air 
after  the  decay  of  the  cloud  so  that  the  convective  heating  and  moistening 
is  proportional  to  the  temperature  and  mixing  ratio  differences  between  the 
cloud  and  the  surrounding  atmosphere. 

The  net  convergence  of  moisture  into  a  vertical  column  of  air  cf  a 
unit  cross  section  produced  by  the  large-scale  flow  and  by  the  evaporation 
from  the  surface  is  expressed  by 


where  g  is  the  gravity,  V  is  the  horizontal  velocity,  q  is  the  water  vapor 
mixing  ratio,  p  is  the  density,  Pg  is  the  surface  pressure,  Cp  is  the  drag 
coefficient  and  subscripts  g  and  0  denote  the  surface  and  anemometer  levels. 
The  amount  of  moisture  convergence  needed  to  create  a  deep  cumulus  cloud  of 
unit  area  and  pressure  depth  ( Pb  -  Pt)  is 


M 

c 


fPb  C 

-  (T  -T) 

g  Jp  L  c 


where  subscript  C  denotes  cloud  variables,  P^  and  P,.  are  the  pressures  at 
the  base  and  top  of  tne  cloud,  L  is  the  latent  heat  of  condensation,  ?  is 
the  specific  heat  at  constant  pressure  and  T  is  the  temperature. 

The  variable  has  two  components:  one  is  the  moisture  required  to 
raise  tne  cloud  temperature  to  T  and  the  other  to  raise  the  cloud  mixing 
’’  t'.i  c  to  q,, .  Therefore, 


where 


and 


C 

t-£  (T  -T)dP 

Li  C 


I 


(qc-q)dP 


(5) 


The  fractional  area  of  a  grid  square  that  will  be  covered  by  newly  formed 
convective  clouds  in  a  time  interval  Ax  is 


a  =  Mt-Ax/Mc  (6) 

The  variable  Ax  can  be  considered  to  be  the  cloud  time-scale  parameter 
TKr ishnamurti  et  al. ,  1983).  In  practice,  the  time  interval  of  a  large- 
seal  ■?  prediction  model,  At,  is  commonly  used  in  place  of  Ax. 

Since  Kuo's  scheme  stipulates  that  the  cloud  will  mix  its  temperature 
and  mixing  ratio  with  those  of  the  environmental  air,  the  large-scale  tem¬ 
perature  T  and  mixing  ratio  q  after  mixing  are  expressed  by 

Tt  +  A t  =  T#  +  a(Tc  "  ^ 
q’  +  a(qc  -  q*) 


I 
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where  T  is  the  temperature  at  time  t  At  before  cumulus  parameterization 
and  q'  is  the  mixing  ratio  at  the  previous  time  step.  In  this  process,  all 
of  the  moisture  convergence  is  used  to  build  clouds,  leaving  the  environ¬ 
mental  mixing  ratio  unaffected  by  the  large-scale  flow.  In  other  words, 

* 

a(qc~q  )  in  Kuo's  scheme  includes  the  terms  of  large-scale  advection,  sur¬ 
face  evaporation,  cloud  condensation  and  cloud  transport  of  moisture.  The 
combination  of  these  processes  will  moisten  the  atmosphere  and  is  commonly 
called  the  cloud  moistening  effect  in  Kuo's  scheme.  When  cloud  condensation 
and  cloud  transport  alone  are  considered,  they  will,  in  general,  dry  the 
atmosphere  and  are  called  the  cloud  drying  effect.  To  compute  the  cloud 

drying  effect  using  Kuo's  scheme,  large-scale  advection  and  surface  evapora- 

* 

tion  should  be  subtracted  from  a(qc-q  )  profiles. 

2.2  Kuo  (197*0 

Kuo's  1965  scheme  of  cumulus  parameter izat ion  has  been  remarkably 
successful  in  hurricane  simulation  studies  (Rosenthal,  1970).  However,  it 
underestimated  the  convective  rainfall  and  heating  rates  in  large-scale 
tropical  applications.  This  defect  is  primarily  related  to  the  allocation 
of  Mt  into  IQ  and  I  .  Since  Iq  is  typically  twice  as  large  as  Ig  in  the 
tropics,  Kuo's  1965  scheme  uses  about  1/3  of  the  moisture  supply  to  heat  the 
atmosphere  and  2/3  to  moisten  the  atmosphere.  Observations  show  that  almost 
all  of  the  moisture  supply  goes  to  heat  the  atmosphere.  Kuo  (197*1)  recog¬ 
nized  this  shortcoming  and  revised  the  scheme  to  provide  a  more  reasonable 
subdivision  of  this  moisture  supply. 

In  deriving  the  new  scheme,  the  equations  for  the  potential  tempera¬ 
ture  9  and  the  water  vapor  mixing  ratio  q  of  the  large-scale  system  are 
written  in  p-coordinates  as 
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d0  _  _ L_  -  _  _L_  _  am'  6' 

dt  r  c  n  c  tt  3p 
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where  o’  and  ~c*  are  the  condensation  rates  produced  by  large-scale  and  cloud 
convective  motions,  respectively,  and  Qp  is  the  heating  rate  by  radiation 
and  turbulent  diffusion.  is  the  rate  of  mixing  ratio  change  by  turbulent 
diffusion  and  it  is  the  non-dimensional  pressure  defined  a3: 


IT 


(P/P  ) 
0 


R/C 

P 


The  terms  on  the  right  hand  sides  of  (7)  and  (8)  are  the  cumulus  effect  on  e 
and  q. 

In  order  to  correct  the  disproportion  of  Mt  into  IQ  and  1^,  Kuo  pro¬ 
posed  that  the  fraction  (1-b)  of  the  moisture  supply  Mt  is  condensed  during 
cumulus  convection  and  is  precipitated  out  as  rain  while  the  remaining  frac¬ 
tion  b  of  Mj.  is  stored  in  the  air  to  increase  the  humidity: 
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where  R  is  the  rainfall,  and  b  is  the  parameter  for  partitioning  the  mois¬ 
ture. 

The  rate  of  the  release  of  latent  heat  by  deep  cumulus  convection  due 
to  condensation  can  now  be  written  as: 


aQc 


L  -* 
-  c 

C  IT 

P 


gO-b) 


L 

C  TT 
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,  V!i 

t  <p6-pt) 


(11  ) 


where  Qc  is  the  latent  heating  effect  by  one  unit  area  of  cloud,  Nk9(p)  is 
the  vertical  distribution  function  of  Qc.  Kuo  proposed  that  NkQ(p)  could  be 
determined  by  the  cloud  -  environment  temperature  differences.  Therefore, 
it  can  be  written  as 
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ke 


(P) 


(9-9) 

c 

<0  -0> 
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which  satisfies  the  condition: 


(12) 
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pt 

The  angle  brackets  represent  the  vertical  averaging  operator  which  is  de¬ 
fined  by: 


()  dP/(P  -P,  ) 
b  t 


(13) 
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In  Kuo's  1965  scheme  (previous  section)  entrainment  was  not  considered 
because  only  deep  clouds  were  included.  The  influence  of  entrainment  on  9C 
can  be  taken  into  consideration  either  by  the  employment  of  a  simple  empiri¬ 
cal  formula  or  by  the  use  of  a  simple  entrainment  model.  Kuo  proposed  the 
following  equation  expressed  in  T  instead  of  6: 

J  Y 

T  -T  =  T  -[l+I  0,(1  -  P/P  )  J]T  (14) 
c  3  j-T  J  0 

where  Ts  is  the  temperature  along  the  moist  adiabatic  process,  the  summation 
is  over  the  different  types  of  cumulus,  Oj  is  a  constant  that  is  determined 
by  the  vertical  extent  of  the  particular  type  of  cumulus,  and  y,  is  another 
constant  which  is  most  probably  2.  Similar  equations  can  be  used  to  com¬ 
pute  the  entrainment  effect  on  the  cloud  mixing  ratio  (qc~q).  Kuo  also 
proposed  a  simple  formula  to  compute  the  cloud  vertical  velocity  which, 
together  with  Tc  and  qp  computed  from  (14),  can  be  used  to  compute  the  ver¬ 
tical  transports  of  heat  and  moisture  by  cumulus  clouds  [see  (7)  and  (8)]. 
Kuo  argued  that  the  vertical  sensible  heat  transport  is  much  smaller  than 
aQc  and  that  its  effect  is  only  to  shift  the  level  of  maximum  heating  rate 
to  a  slightly  higher  altitude.  The  computation  of  c"*  and  the  flux  diver¬ 
gence  of  heat  and  moisture  complete  Kuo's  1974  cumulus  parameterization. 

In  the  above  formulations  Kuo  did  not  give  a  formal  expression  for  the 
b  parameter  to  close  the  system.  He  only  assumed  that  b  is  much  smaller 
than  1  in  regions  of  low-level  convergence  in  the  tropics,  which  implies 
that  all  of  the  moisture  convergence  brought  about  by  the  large-scale  flow 
is  precipitated  out  as  rain. 
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2.3  Anthes  (1977) 

Following  Kuo’s  197^  scheme,  Anthes  (1977)  argued  that  the  vertical 
distribution  of  heating  by  condensation  is  not  adequately  represented  by  the 
cloud-environment  temperature  differences.  Therefore,  he  proposed  to  use 
the  actual  cloud  condensation  profile  from  a  cloud  model  for  the  vertical 
heating  distribution  function: 

NAe(P)  •  C*/<C*>  (15) 

The  vertical  distribution  of  condensational  heating  is  therefore  given  by 


aQc 


_k_  8( 1  M  n  /p) 

c  IT  (P  -  P„)  t  A9  “ 
p  b  t 


(16) 


which  is  the  same  as  (11)  except  for  the  vertical  heating  distribution  func¬ 
tion. 

He  also  proposed  that  b  is  a  function  of  the  mean  relative  humidity  in 
the  column,  increasing  from  near  zero  in  a  saturated  column  to  1  as  the  mean 
relative  humidity  approaches  a  critical  value  RHC>  Therefore  b  is  given  by 


b 


1 -<RH> 
1-RH 


<RH>  >  RH 


1  <RH>  <  RH 

c 


(17) 


where  n  is  a  positive  exponent  of  order  1  which  may  be  empirically  adjusted. 
The  moisture  equation  at  gridpoints  with  convection  is  written  as 


where  3q  /3t  is  the  combined  moistening  effect  due  to  water  vapor  convei — 


genee  by  the  large-scale  flow  and  the  removal  of  water  vapor  by  condensa¬ 
tion.  Given  a  value  of  b,  the  vertical  integration  of  3q/3t  is 
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bM 
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He  proposed  the  following  form  for  the  vertical  distribution  of  3q+/3t: 


H-  ■  "  «  a5P)  (,9> 

where  Mftq(p)  is  the  vertical  moistening  distribution  function.  Concluding 
that  the  most  moistening  should  occur  in  the  driest  layers,  he  suggested  the 
following  form  for  Nftq(p) 


( 1 00%  -  RH)qs(T) 
NAq(P)  =  <(100$  -  RH)  qg(T)> 


(20) 


where  qg(T)  is  the  saturation  mixing  ratio  at  the  environmental  temperature 
T. 

In  order  to  close  the  cumulus  parameter! zation  scheme  by  using  the 

above  formulation,  a  cloud  model  is  required  to  provide  the  mean  cloud  prop- 
„  _  * . 

erties  Vo)c,  T  ,  qc  and  c  )  for  computing  the  cloud  transport  terms  and  the 
vertical  heating  distribution  function  NA0(p).  In  this  scheme,  the  vertical 
distribution  of  the  latent  heating  is  very  sensitive  to  the  radius  of  the 
mean  cloud.  How  to  determine  the  radius  of  the  mean  cloud  is  the  major 


problem. 


2.4  Krlahnamurtl  et  al.  (1976) 


Krishnamurti  et  al.  modified  Kuo's  1974  scheme  and 
hypothesis  for  computing  b.  First,  they  used  the  following 
and  moisture  equations: 
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(22: 


They  defined  the  available  moisture  supply  as  the  vertically  integrated  ve 
tical  advection  of  mixing  ratio: 


I  f  6 

g  Jn 


w  fp  dP 


(23) 


Furthermore,  they  defined  Qq  and  Q„  as 

o  q 


a  0 
0  0 


‘8  f  b  °p  T  ,V 


de- 


0  I  _P  T  r  . 

T  JD  l  e  l—  +  “  apJ  dp  "  (1‘b)I 


(24) 


a  Pb 

Vq  =  f  Jp  (Vq)/AT  dP 


bl 


(25) 


+  Q 

w0  wq 


where  Q  is  a  measure  of  the  total  moisture  supply  needed  to  cover  a  unit 
horizontal  area  with  clouds,  is  the  moisture  supply  required  to  raise  the 
temperature  to  the  cloud  temperature  in  this  column  and  0  is  the  moisture 
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supply  required  to  saturate  the  column  at  the  cloud  temperature.  According 
to  their  argument,  the  inclusion  of  the  second  term  on  the  righthand  side  of 
(24)  provides  a  smooth  transition  between  the  convective  heating  in  the 
conditionally  unstable  atmosphere  and  the  large-scale  condensational  heating 
in  the  stable  situation.  The  detailed  description  of  the  smooth  transition 
when  this  term  is  retained  can  be  seen  in  Krishnamurti  et  al.  (1983).  The 
cloud  temperature  and  mixing  ratio  are  computed  from  a  moist  adiabatic  proc¬ 
ess  so  ec  and  qc  in  (24)  and  (25)  assume  their  saturated  moi3t  adiabatic 
values  93  and  q  ,  respectively.  The  following  hypothesis  is  made  to  close 
the  system: 


1  f 
g  J 


J?  T  f 30  + 

l  e  k  at 


v-ve)  dp 


(f*  *  ^'Vq)  dP 


(26) 


The  implication  here  is  that  in  the  cloud  column  the  total  change  of  9  and  q 
due  to  both  vertical  advection  and  convection  are  proportional.  This  condi¬ 
tion  brings  the  9  and  q  values  to  the  moist  adiabatic  values  in  near  unison. 
Substituting  (21)  and  (22)  into  (26)  and  using  A  defined  as 
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we  can  get 
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Once  b  is  determined,  aQ  and  aq  can  be  computed  from  (24)  and  (25).  These 
values  in  turn  can  be  used  to  compute  the  cumulus  effect  in  (21)  and  (22). 
In  notation  similar  to  Kuo  (1974),  the  convective  heating  profile  can  be 
expressed  as 


aQc 


L  gd-tOH.  "  KRkP> 


C  7T 
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where  the  vertical  heating  distribution  function  is 
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The  mixing  ratio  equation  can  be  written  a3 
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Basically  the  above  approach  is  the  extension  of  the  Kuo  (1974)  scheme.  The 
closure  assumption  is  ad  hoc  and  is  not  based  on  any  valid  theory. 

2.5  Molinarl  (1982) 

Based  on  the  framework  developed  by  Krishnamurti  et  al.  (  1976),  an 
alternative  closure  was  proposed  by  Molinari  (1982)  for  the  b  parameter 
which  more  explicitly  requires  the  warming  and  moistening  of  the  cumulus 
column  to  produce  a  smooth  change  toward  moist  neutrality  in  both  tempera¬ 
ture  and  mixing  ratio.  The  closure  is  written  as 


t>  f  b  f  b  r  b 
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Substituting  (21)  and  (22)  into  (32), 


yields 
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(33) 


This  procedure  forces  the  instantaneous  time  tendencies  of  0  and  q  to  be 
such  that  they  would  reach  their  moist  neutrality  simultaneously. 

This  scheme  is  an  extension  of  that  of  Krishnamurti  et  al .  (1976).  The 
only  difference  between  these  two  schemes  is  the  closure  assumption.  Accord¬ 
ing  to  Molinari's  (1982)  results,  the  proposed  closure  produces  a  more  accu¬ 
rate  time  variation  of  stability  that  approaches  moi3t  neutrality  better 
than  Krishnamurti  et  al.  (1976).  However,  the  closure  assumption  for  b  is 
still  an  ad  hoc  assumption. 

2.6  Krishnamurti  et  al.  (1983) 

Following  the  previous  work  of  Krishnamurti  et  al.  (1976),  they  intro¬ 
duced  an  additional  source  of  moisture  supply,  In,  which  may  be  considered  a 
nonmeasurable  mesoscale  (or  subgrid-scale)  supply.  The  total  supply  of 
moisture  is  expressed  by 


I '  *  I ( 1 +n) 


(3*0 


and  the  rainfall  rate  R  and  moistening  rate  M  are  given  by 


-15- 


R  =  I’(l-b)  =■  I ( 1 +n) ( 1 -b)  (35) 

M  »  I  ’  b  I( 1 +n)b  (36) 

In  this  formulation  there  are  two  unknown  parameters,  n  and  b,  that  need  to 
be  determined  for  the  closure  of  the  cumulus  parameterization  scheme.  They 
proposed  that  n  and  b  be  expressed  as  functions  of  pairs  of  known 

large-scale  variables  such  as  the  vertically  averaged  large-scale  vertical 
velocity  w  and  the  700  mb  relative  vorticity  £  .  Based  on  the  multiple  re¬ 
gression  analysis,  they  get 
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(90) 


They  have  tested  the  above  scheme  to  determine  the  vertical  distribu¬ 
tions  of  the  heating,  moistening,  and  rainfall  rates  for  the  entire  GATE 
('Global  Atmospheric  Research  Program’s  Atlantic  Tropical  Experiment)  Phase 
III  data.  The  results  of  these  calculations  are  in  good  agreement  with  the 
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observed  estimates.  However,  they  did  not  mention  whether  the  regression 
coefficients  are  universal  constants  or  not.  Those  constants  may  change 
with  different  time  periods  or  locations. 

3.  Semi  prognostic  Teat  of  the  Kuo  Scheme 

Imperfectly  represented  cloud  processes  in  cumulus  parameterization 
schemes  may  interact  with  the  large-scale  flow  in  unrealistic  ways.  For 
this  reason,  any  proposed  parameterization  scheme  must  be  verified  with  ob¬ 
served  data.  Generally  speaking,  two  approaches  have  been  taken  to  test  a 
cumulus  parameter i zat ion  scheme  against  observations.  The  first  approach 
may  be  categorized  as  a  semi  prognostic  approach  because  it  is  a  one  time- 
step  prediction  of  cumulus-produced  heating  and  drying  profiles  and  rainfall 
rate.  The  other  approach  is  to  incorporate  the  parameterization  schemes 
into  a  large-scale  prediction  model,  such  as  a  regional,  hemispherical  or 
global  model. 

In  this  section  we  will  focus  on  a  semi  prognostic  test  of  the  Kuo 
scheme  using  GATS  B-scale  data.  The  apparent  heat  source  Q1  and  the  appar¬ 
ent  moisture  si.nK  Q?  are  defined  by  the  following  equations: 


1  rli 

it  "  at 


V  •  70 


39 1 
w  3P ' 


l  r3ui'e'  I 

IT  k  3P  J 


) 


3 


2 


V  •  7q 


w 


111 

3P  J 


3u)’q’  -I 

3P  J 


(42) 


Here  the  overbar  denotes  the  average  over  the  GATE  B-scale  area  and  tne 
prime  denotes  the  cloud  scale.  Q,-QR  represents  the  cloud  heating  effect 
and  3 ->  represents  the  cloud  drying  effect.  Dr.  Y.-L.  Chen  of  the  University 
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of  Hawaii  has  diagnostically  estimated  Q1  and  Q?  for  the  entire  GATS  Phase 
III  (30  August  -  IS  September)  B-scale  area  by  computing  the  terms  between 
the  equal  signs  in  (41)  and  (42).  We  will  use  his  data  and  refer  to  his  Q, 
and  Qp  values  as  the  observation  values.  The  vertical  resolution  of  the 
data  is  33  uneven  levels  from  the  surface  to  92.5  mb.  The  area  average  is 
over  2.5°  in  radius  and  the  time  resolution  is  every  3  hours.  The  radi- 
ational  heating  (or  cooling)  rates  QR  that  we  used  in  this  study  are  taken 
from  Krishnamurti  et  al.  ( 1  9 8 3 ) .  Table  1  lists  the  6-hourly  net  radiational 
heating  rates  averaged  over  the  entire  phase  III  period.  In  this  data, 
significant  changes  in  the  magnitude  of  QR  are  evident  between  the  night  and 
daytime  hours. 

In  the  s“m i prognost i c  approach,  the  cumulus  heating  and  drying  ef¬ 
fects,  i.e.,  those  terms  on  the  righthand  side  of  (41)  and  (42)  excluding 
^P,  are  calculated  from  the  parameterization  schemes  (thereafter  referred  to 
is  the  calculated  values)  and  compared  with  the  observation  values  of  Q 1 - 
and  -i?.  The  ,  Ivan*,  age  of  this  approach  is  that  a  prediction  can  be  made 
without  using  a  large-scale  prediction  model.  Moreover ,  since  this  is  only  a 
one  time-step  int»gr  it  ion,  the  predicted  results  are  free  from  numerical 
nonlinear  errors  md  complicated  physical  feedback  effects.  Any  successful 
p  ir  meter  i  cat  i  s 'hem-1  ha  .a  to  first  achieve  a  good  score  on  this  semi- 
P"  '  gr  os  t  i  1  . 

Ir.e  method  used  to  parameter  i  ze  the  cumulus  effects  in  this  test  is  a 
’-mini- ;c,  of  K  <  o  ’  .3  ’  965  ‘.rid  '974  schemes.  They  can  be  summarized  as  foi- 

>  •  In*-  m-  ripply  is  .wr.puf  **d  from  t n*’  vertically  integrated 

m-  :  "  ,’•*  •.  .-r md  surf  ice  ev  tporatior. . 
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Table  1 .  Six-hourly  total  radiative  heating  rate  (C/day)  averaged  over  the 
entire  Phase  III  of  the  GATE  period. 


Pressure 
level (mb ) 

0000-0600 

0600-1200 

1200-1800 

1800-2400 

100 

-0.7S 

-0.09 

-0.10 

-0.74 

1  ”  5 

-0.95 

-0.14 

-0.15 

-0.92 

250 

-1.52 

-0.51 

-0.54 

-1.55 

525 

-2.02 

-0.29 

-0.56 

-2.14 

400 

-2.54 

-0.40 

-0.47 

-2.45 

4  7  5 

-2.55 

-0.53 

-0.S9 

-2.52 

550 

-2.70 

-0.76 

-0.70 

-2.56 

625 

-2.56 

-0.79 

-0.72 

-2.24 

700 

-  2.  OS 

-0.89 

-0.35 

-  2.01 

775 

-1.31 

-0.95 

-0.39 

-1.79 

850 

-1.46 

-0.80 

-0.79 

-1.46 

925 

-1.04 

-0.49 

-0.48 

-1.05 

1000 

-0.10 

-0.10 

-0.12 

-0.12 

.w.v-/ 


(b)  A  fraction  (1-b)  of  the  moisture  supply  is  used  to  heat  the  at¬ 
mosphere  and  a  fraction  b  of  the  moisture  supply  is  used  to  mois¬ 
ten  the  atmosphere. 

(c)  The  calculated  (Q1~QB)(total  cloud  heating)  profile  is  computed 
using  (11). 

(d)  The  cloud  moistening  effect  (including  the  large-scale  advection 


and  surface  evaporation)  is  computed  from 


3a  8b*t  (V<1) 

3t  '  <W  <V<1> 


(e)  The  calculated  Q?  profile  is  computed  by  subtracting  the  large- 
scale  advection  and  surface  evaporating  terms  from  the  above 
equation. 

The  partition  of  the  moisture  supply  is  based  on  Kuo's  1 97 ^  scheme.  Once 
this  partition  is  determined,  the  Q^Qr  and  Q2  profiles  are  computed  based 
on  Kuo ' 3  1965  scheme. 

As  mentioned  in  Section  2.2,  Kuo  proposed  that  the  parameter  b  in  his 
1974  scheme  is  very  small.  Several  methods  to  compute  b  have  been  proposed 
(see  Section  2)  but  no  verification  has  been  made.  However,  Krishnamurti  et 
al.  (1980)  compared  the  observed  rainfall  to  that  predicted  from  Kuo's 
scheme  with  small  b  (b  *  0)  and  found  extremely  good  agreement.  Therefore, 
in  this  semiprognostic  test  we  have  set  b  equal  to  zero.  For  simplicity, 
the  entrainment  effect  on  cloud  temperature  and  mixing  ratio  is  not  in¬ 


cluded. 
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Figs.  1  and  2  are  the  calculated  and  observed  Q1  -  QR  for  September 
1-18,  197*1  shown  on  pressure-time  sections.  The  large  observed  heating 
occurred  on  the  following  dates:  2,  *4,  5,  12,  1 3 >  1**.  15,  16,  and  17.  Tnere 
is  general  agreement  between  the  calculated  and  the  observed  values.  To  ex¬ 
amine  the  results  closely,  Fig.  3  shows  the  comparison  of  the  calculated 
Q1  -  Qr  profiles  with  those  obtained  from  observations  during  periods  of 
heavy  rainfall  on  2(1200,  1800GMT),  4(1200,  1800GMT),  5(1200,  1800GMT), 
12(1  200GMT ) ,  1 5 ( 0000GMT ) ,  16(1200GMT),  and  17(1200GMT).  The  maximum  value 
of  heating  lies  between  10  and  15°C/day,  in  good  agreement  with  the  observed 
heating  values. 

Fig.  4  shows  the  comparison  of  calculated  and  observed  Q1  -  QR  during 
periods  of  weak  convection  period  on  2(0600GMT),  3(1800GMT),  9(0000GMT), 
10(1 300GMT ) ,  11(1 200GMT) ,  12(0000GMT),  13O200GMT),  14(1800GMT),  and 
16(0600,  1800GMT).  The  observed  maximum  heating  is  located  below  550  mb; 
however,  the  calculated  maximum  heating  is  above  550  mb.  The  cloud  heating 
rate  has  been  consistently  overestimated  above  500  mb  and  underestimated  at 
the  lower  levels. 

In  general,  the  observed  vertical  distributions  of  convective  heating 
are  highly  variable;  but  the  vertical  profiles  of  the  cloud-environment 
temperature  difference  computed  from  rising,  undiluted  surface  parcels  (Tc  - 
T)  are  less  variable  (Song  and  Frank,  1983).  Since  the  cloud  heating  pro¬ 
file  was  based  on  the  profile  of  Tc  -  T  and  no  entrainment  was  considered, 
the  poor  agreement  between  the  calculated  and  observed  cloud  heating  pro¬ 
files  during  weak  convective  conditions  is  not  surprising.  On  the  other 
hand,  during  strong  convection  the  cloud  heating  profile  does  match  the 
undiluted  Tp  -  T  profile  well. 
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The  time-averaged  (for  the  entire  GATE  Phase  III)  vertical  profiles  of 
Q^-Qp  are  shown  in  Fig.  5.  Again,  the  calculated  Q,  -  QR  was  overestimated 
above  700  mb  and  underestimated  under  700  mb.  The  maximum  difference  be¬ 
tween  them  is  about  2.5°C/day. 

4.  Present  Approach 

The  results  from  the  above  semi  prognostic  tests  indicate  that  Kuo's 
scheme  would  produce  good  results  during  periods  of  deep  clouds  and  heavy 
precipitation.  However,  it  also  indicates  the  need  to  include  the  entrain¬ 
ment  according  to  convective  activity  in  computing  the  cloud  heating  pro¬ 
file.  The  entrainment  should  have  a  minimal  effect  on  the  profiles  of  Tc  -T 
during  strong  convection;  but  it  will  significantly  alter  the  TQ  -  T  pro¬ 
files  and  lower  the  level  of  maximum  cloud  heating  during  weak  convection, 
as  indicated  by  observations. 

When  Kuo  proposed  his  1974  scheme,  he  actually  formulated  the  effect 
of  entrainment  on  cloud  heating  profiles  [see  (14)].  This  part  of  his 
scheme  has  received  little  attention  since  it  requires  that  two  additional 
constants  be  specified.  Otherwise,  his  entrainment  consideration  is 
straightforward  and  easy  to  apply.  Our  improvement  of  Kuo’s  scheme  includes 
the  adoption  of  the  entrainment  formulation  based  on  Kuo's  method,  the  de¬ 
tails  of  which  will  be  discussed  later. 

Other  than  entrainment,  Kuo's  1974  scheme  requires  the  computation  of 
the  parameter  b,  the  fraction  of  moisture  convergence  used  to  moisten  the 
atmosphere.  The  value  of  1-b  is  the  fraction  of  moisture  convergence  pro¬ 
ducing  precipitation.  Sui  and  Yanai  (1984)  found  a  remarkable  similarity 
between  the  time  variations  of  the  percentage  coverage  of  deep  clouds  with 


tops  above  the  300  mb  level  and  the  rainfall  rate.  The  rainfall  rate  can 
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therefore  be  considered  proportional  to  the  activity  of  deep  cumulus  clouds. 
Chao  et  al .  (’979)  also  pointed  out  that  the  cloud  tops  were  substantially 
higher  during  the  disturbed  period  than  during  the  undisturbed  period.  Song 
and  Frank  (1983)  suggested  "an  area  for  future  research  would  be  to  detei — 
mine  whether  a  relationship  exists  between  the  area-a veraged  rainfall  rate 
and  mean  cloud  entrainment  rate  which  would  bring  profiles  of  TQ  -  T  and  net 
convective  heating  into  closer  agreement".  These  papers  indicated  that 
(l-b)  could  be  used  as  an  indicator  of  cloud  activity  and  the  entrainment 
effect  on  the  model  cloud.  We  have  used  the  GATE  Phase  III  data  to  calculate 
the  b  value  according  to  (28)  and  (33).  The  results  indicate  that  the  b 
value  from  these  two  equations  are  almost  identical  and  that  their  values 
are  very  reasonable  in  comparison  with  precipitation  observations.  There¬ 
fore,  (28)  will  be  used  to  determine  the  b  value. 

In  order  to  use  (19),  as  Kuo  proposed  for  considering  the  influence  of 
entrainment,  the  value  of  aj  and  Yj  has  to  be  predetermined .  At  this  mo¬ 
ment,  we  assume  that  only  one  particular  type  of  cumulus  cloud  exists  at 
each  time  step.  Therefore,  only  a1  and  Y1  need  to  be  determined  before  (19) 
is  U3ed.  Since,  according  to  Kuo,  the  value  of  Y1  is  most  likely  greater 
than  2,  the  value  of  Y1  in  this  study  is  set  empirically  to  2.2.  Because  a1 
is  U3ed  to  determine  the  vertical  extent  of  the  cumulus  cloud,  we  need  to 
define  the  level  of  cloud  top,  Pj,  where  the  cloud  heating  equals  zero.  From 
several  tests,  we  found  that  the  best  equation  to  estimate  the  vertical 
extent  of  the  cumulus  cloud  is 


P 


T 


900  x  EXP  (-1 .5(1-b)2} 


(93) 
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Fcr  periods  of  strong  convection  (b  •>  0),  Pj  *  200  mb  while  for  periods  of 
extremely  weak  convection  (b  *■  1),  Pj  *  900  mb.  Once  P-j.  is  determined  for 
the  particular  cumulus  cloud  under  consideration,  a-]  is  computed  using  (14) 
by  setting  Te  -  T  =  0  at  the  cloud  top,  Pp  This  computation  leads  to 


Let  us  consider  two  cases,  each  with  the  same  Mt  and  Nkg(p).  In  one 
case,  the  atmosphere  is  stable  and  Tc  -  T  is  small.  In  the  other  case,  the 
atmosphere  is  unstable  and  Tc  -  T  is  large.  From  (11)  and  (12),  the  two 
cases  should  produce  the  same  cloud  heating  profile,  which  is  not  realistic 
since  the  stability  is  not  considered  in  the  determination  of  the  cumulus 
activity.  Song  and  Frank  (1983)  found  that  the  net  conditional  instability 
(Cl)  seldom  dropped  below  1 °C  during  the  convection,  where  Cl  was  defined  as 


For  this  reason,  we  set  a  critical  Cl  value  (Cl  cr )  as  1  °C .  If  Cl  >  CInr, 
there  is  no  cloud  heating.  The  entire  procedure  of  this  cumulus  parameter¬ 
ization  scheme  is  rresented  in  Fig.  6  in  flow-chart  form. 
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Present  Scheme 


Seaiprognostic  Test  of  the  Present  Scheme 


Using  the  present  cumulus  parameterization  scheme,  we  again  calculated 
Q-)  ~  QR  profiles  for  the  entire  GATE  Phase  III  period.  During  periods  of 
strong  convection,  the  results  are  no  different  than  the  previous  computa¬ 
tion.  This  can  be  expected  since  the  entrainment  effect  is  small  for  deep 
convection.  Fig.  7  shows  the  cloud  heating  profiles  using  the  present 
scheme  and  observed  Q1  -  QR  during  selected  periods  of  moderate  to  weak  con¬ 
vection,  as  in  the  last  section.  We  can  see  a  significant  improvement  com¬ 
pared  to  Fig.  4.  The  level  of  maximum  heating  has  shifted  from  the  upper 
troposphere  to  a  lower  level.  The  cloud  tops  vary  with  convective  activity 
and  can  be  as  low  as  500  mb.  In  general,  the  calculated  maximum  heating  in 
the  present  scheme  is  about  the  same  as  in  the  observations. 

Fig.  3  shows  the  pressure-time  sections  of  the  calculated  Q1  -  QR 
fi  Id  for  September  1-18,  1 97 ^  in  the  GATE  Phase  III  period.  A  comparison 
of  Fig.  8  with  Fig.  1  shows  the  following  features:  on  6  September  the 
maximum  heating  has  increased  to  6°C/day  and  was  around  575  mb;  on  12  Sep¬ 
tember  the  maximum  heating  shifted  from  above  400  mb  to  below  500  mb;  and  on 
16-17  September  the  maximum  heating  was  around  500  mb  instead  of  400  mb,  a3 
in  Fig.  1 .  A  comparison  of  Fig.  8  with  Fig.  2  shows  that  the  results  from 
the  present  scheme  agree  well  with  the  observations.  Overall,  the  heating 
distribution  from  the  present  scheme  is  quite  impressive. 

Figs.  9  and  10  show  the  pressure-time  sections  of  the  calculated  and 
observed  apparent  moisture  sink  Q2  field.  The  observed  field  of  Q2  exhibits 


large  drying  >  0)  on  the  following  dates,  2,  4,  5,  12,  13,  14,  16,  and 

17.  The  maximum  intensity  of  this  drying  occurs  around  and  below  700  mb. 
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Fig.  9.  Vertical  cross  section  of  the  calculated  cloud  drying  effect  ( Q2 ) 
for  the  present  scheme. 
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Occasionally ,  it  extends  to  very  high  levels.  The  periods  of  strong  drying 
usually  coincide  with  the  periods  of  strong  cloud  heating.  In  general,  the 
present  scheme  essentially  captures  all  of  the  observed  features  of  Q2- 

The  time  averaged  vertical  profiles  of  calculated  Q1  -  and  Q2  are 
very  similar  to  the  observed  profiles  (see  Figs.  11  and  12).  The  maximum 
heating  level  is  around  600  mb  in  both  profiles,  with  the  maximum  difference 
between  them  being  less  than  1.1°C/day.  The  calculated  Q2  profile  is  almost 
coincident  with  the  observed  profile.  For  a  comparison,  Fig.  13  shows  the 
same  vertical  profiles  of  Q1  -  and  Q2  as  calculated  by  Krishnamurti  et 
al.  (1983). 

Fig.  14  is  a  time  series  of  rainfall  rates  (mm/day)  from  1-18  Septem¬ 
ber  1974  as  estimated  by  the  present  scheme,  as  well  as  those  from  the 
moisture  budget  and  the  radar  measurements.  The  estimated  rainfall  from  the 
moisture  budget  Rbu  13  calculated  from: 


e 

L  g  g 


where  E  is  the  estimated  surface  evaporation  rate  calculated  using  the  bulk 

o 

aerodynamic  formula: 


Eg  -  OoVo'VV  ("8) 

where  q  is  the  surface  mixing  ratio  calculated  based  on  the  sea  surface 

o 

temperature  given  by  Wuchnitz  et  al.  (1977).  The  calculated  rainfall  rate 
agrees  closely  with  the  two  independent  estimates  from  the  moisture  budget 
and  from  radar  measurements.  The  peak  calculated  rainfall  values  on  2,  14, 
and  16  September  are  almost  identical  to  the  observed  values.  There  are  some 


Time-averaged  vertical  distribution  of  the  cloud  heating  profiles 
( Q1 — Qr ) -  The  observed  value  (dashed  line)  is  compared  with  those 
calculated  from  the  present  scheme  (solid  line). 
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TIME  AVERAGED  Q2  DURING  PHASE  3  GATE 
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Fig.  12.  Time-averaged  vertical  distribution  of  the  cloud  drying  profiles 
(Q2).  The  observed  valued  (dashed  line)  is  compared  with  those 
calculated  from  the  present  scheme  (solid  line). 
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rig.  1U.  A  time  series  of  rainfall  rates  from  1-13  Septemoer  estimates  Py 
the  present  scheme,  the  moisture  budget  and  the  radar  measure¬ 


ments  . 


differences  in  the  time  of  maximum  calculated  rainfall  rate  and  the  observed 
maximum  by  the  radar  during  the  periods  of  5-6  and  9*10  September.  No  con¬ 
clusive  explanation  can  be  made  with  the  limited  amount  of  data. 

In  order  to  examine  the  general  behavior  of  the  present  scheme  in 
various  conditions  in  the  tropics,  the  composite  easterly  wave  data  in  Phase 
III  were  also  used  in  this  semiprognostic  study.  The  observed  Q1  values 
were  provided  by  Dr.  Chen  of  the  University  of  Hawaii.  Fig.  15  shows  the 
calculated  and  observed  Q-,  -  QR  fields  for  the  composite  wave  in  the  GATE 
B-scale  area  at  latitude  8.5°N.  The  waves  are  separated  into  8  categories 
with  category  4  representing  the  trough  area  and  category  8  representing  the 
ridge  area.  We  can  see  that  the  level  of  the  maximum  heating  rate  varies 
with  the  wave  category.  The  comparisons  between  observed  and  calculated  Q1 
-  curves  are  reasonably  good. 

6.  An  Evaluation  of  Kuo's  Cumulus  Parameterization  Scheme  Using  a  Two- 

Dimensional  Cloud  Cluster  Model 

As  stated  previously,  it  is  necessary  for  a  good  cumulus  parameteri za- 
t ion  scheme  to  produce  cloud  heating  and  drying  effects  and  precipitation 
rates  in  agreement  with  diagnostic  results  in  a  semi prognost i c  test.  How¬ 
ever,  this  agreement  does  not  imply  that  these  schemes  will  have  prognostic 
value.  This  point  is  clear  since  a  'scheme'  requiring  no  change  in  the 
large-scale  temperature  and  moisture  fields  will  give  good  verification  in 
the  semi -prognost i c  approach  but  cannot  be  used  in  a  large-scale  prediction 
model . 

In  this  step  of  tne  work,  a  cloud  cluster  model  will  be  used  to  evalu¬ 
ate  the  cumulus  parameterization  scheme  in  a  fully  prognostic  sense.  The 
cloud  cluster  model  with  a  fine  grid  has  been  used  to  simulate  the  develop¬ 
ment  of  clouds  and  a  cloud  cluster  under  a  weak  low-level  lifting  condition 
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Fig.  15.  The  vertical  distributions  of  the  cloud  heating  profiles  (Q.|- 
for  the  composite  waves  during  GATE  Phase  111  period.  Tine  calc 
lated  values  are  indicated  by  solid  lines  and  the  observed  valu 
are  indicated  by  dashed  lines. 
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in  the  GATE  environment  (Soong  and  Ogura,  1982).  The  model  also  provides 
the  cloud  heating  and  drying  effects  in  the  cloud  cluster  region  as  well  as 
the  large-scale  lifting  generated  by  the  cloud  effects  at  each  integration 
time  step.  This  model  is  modified  to  use  a  coarse  grid  and  the  hydrostatic 
system  of  equations.  The  cloud  dynamics  and  physical  processes  are  elimi¬ 
nated  and  the  present  scheme  of  cumulus  parameterization  is  incorporated 
under  the  same  weak  low-level  lifting  condition. 

6.1.  Model  equations 

The  momentum  equations  of  this  two-dimensional  hydrostatic  cloud  clus¬ 
ter  model  in  the  x-z  slab  symmetric  domain  are: 


3u  3u  3u  „  3it 
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where  0  is  the  potential  temperature,  f  is  the  Coriolis  parameter,  Ug  is 
the  x-component  geostrophic  wind,  and  it  is  the  nondimensional  pressure, 
which  is  defined  as: 


R/C 


it  =  (P/P  ) 
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(51  ) 


where  P  is  pressure,  PQ  is  the  reference  pressure  of  1000  mb,  R  is  the 
gas  constant  of  dry  air  and  Cp  is  the  specific  heat  of  dry  air  at  constant 
pressure . 

The  hydrostatic  equation  is: 

3tt 
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(52) 


where  g  is  gravity  and  9  is  the  virtual  potential  temperature  defined  by 
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The  thermodynamic  and  moisture  equations  are: 


i£  ,  .  UM  -  wli  ♦  a 

at  ax  az  c 


-|f  =  -u-j^-w-^+P  (55) 

3t  3x  3z  r 

where  Qc  is  the  adiabatic  heating  rate  due  to  clouds  and  Pp  is  the  produc¬ 
tion  rate  of  water  vapor  due  to  clouds.  The  radiation  effects  are  neglected 
in  the  thermodynamic  equation. 

The  mass  continuity  equation  is  expressed  in  the  anelastic  form: 

iu  +  llfiw  .  o  ,  (56) 

3x  p  3z 

where  p  is  the  density  and  is  a  function  of  height  only. 

Now  consider  a  variable  A(A  =  it,  0,  or  q)  separated  into  the  horizon¬ 
tal  mean  and  the  departure  as: 

A ( x  ,  z  ,  t )  «T(z,t)+A'(x,z,t).  (57) 

The  horizontal  average  variables  represent  the  large-scale  conditions  and 
are  functions  of  time  and  height  only.  The  primed  variables  define  the 
cloud-scale  system.  Substituting  (57)  into  (i*9),  (50),  (52),  (5*0,  and  (55) 
with  the  use  of  (56),  one  can  obtain  a  set  of  governing  equations  in  flux 


form  as  follows: 
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In  deriving  (60)  it  was  assumed  that  the  horizontal  means  satisfy  the  hydro¬ 
static  relation: 


y-  3  it 
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and  ev  was  approximated  by  1TV  .  In  the  derivation  of  (58),  the  9  in  front 
of  the  pressure  gradient  term  was  approximated  by  0"  .  The  variable  wQ  is 
the  large-scale  low-level  lifting  imposed  on  the  model  based  on  observation. 
As  discussed  by  Ogura  et  al.  (1979)  and  Soong  and  Tao  (1980),  large-scale 
low-level  lifting  usually  preceeds  the  development  of  deep  clouds.  It  will 
destabilize  the  atmosphere  and  provide  the  moisture  supply  for  clouds. 

6.2.  Numerical  technique,  boundary  and  Initial  conditions 

As  shown  in  Fig.  16,  this  model  has  a  horizontal  domain  of  1280  km, 
with  600  km  of  grids  at  the  center  part  of  the  domain  (cloud  cluster  area). 
On  each  side  of  the  cloud  cluster  area,  there  is  an  environment  area.  The 
horizontal  domain  has  64  uniform  gridpoints  with  a  grid  distance  equal  to  20 
km.  It  has  16  grid  intervals  in  the  vertical  with  an  even  grid  spacing  of 
900  m  such  that  the  top  of  the  domain  is  at  the  level  of  14.4  km,  which  is 
around  1 4 0  mb  in  the  P-coordinate . 

A  staggered  grid  arrangement  [the  same  one  used  by  Soong  and  Ogura 
(1973)],  shown  in  Fig.  17,  is  used  in  this  model.  A  leapfrog  time  integra¬ 
tion  and  a  second-order-centered  space  difference  scheme  is  applied  with  a 
time  step  of  60  s.  A  time-smoothing  scheme  formulated  by  Robert  (1966)  is 
adopted  in  the  integration  with  a  smoothing  coefficient  of  0.1  to  avoid  the 
time  splitting.  Thi3  model  also  includes  a  horizontal  second-order  numeri¬ 
cal  smoother  which  is  applied  to  all  the  prognostic  fields.  This  horizontal 
smoother  reduces  the  growth  of  nonlinear  instabilities  and  filters  out  very 
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At  the  upper  and  lower  boundaries,  a  free-slip  boundary  condition  is 
used  for  u  and  a  fixed  boundary  condition  with  zero  velocity  is  used  for 
w.  No  fluxes  of  0  and  q  are  allowed  across  these  boundaries.  The  lat¬ 
eral  boundaries  are  open  (see  Klemp  and  Wilhelmson,  1978).  The  values  of  u 
and  v  at  the  lateral  boundaries  are  determined  by  advecting  corresponding 
values  at  the  adjacent  interior  gridpoints  outward  with  the  gravity  wave 
speed.  As  for  0  and  q  ,  zero  horizontal  gradient  is  assumed  at  the  lat¬ 
eral  boundaries.  The  boundary  conditions  can  be  summarized  as 


w  =  0 


at  upper  and  lower  boundaries 


and 

3u  __  _  _3u 

at  '  cg  ax 

=  -  c  —  at  left  and  right  boundaries  (64) 

at  g  ax 

1!  ,  la  =  o 

3x  3x 

where  C  represents  the  outward-propagating  gravity  wave  speed.  Under  the 

O 

GATE  environment,  the  estimated  gravity  wave  speed  is  36  m  s  1  and  this 
value  is  used  for  C  . 

O 

At  time  t  =  0,  T  and  q  are  horizontally  uniform  and  their  profiles 

are  shown  in  Fig.  18.  These  are  GATE  soundings  taken  from  the  ship  Re¬ 

searcher  at  00  GMT  12  August  1974.  A  computed  low-level  vertical  velocity 
at  00  GMT  12  August  1 97 4  within  an  ITCZ  rainband  (Ogura  et  al.,  1979)  using 

GATE  A/B  scale  data  is  given  in  Fig.  19.  It  is  seen  that  the  maximum  magni¬ 
tude  is  about  1.7  cm  s  1  around  850  mb.  A  sinusoidal  function  was  used  witn 

values  of  3.4  cm  s_1  at  the  center  of  the  domain  and  zero  at  the  boundaries 

of  tne  cloud  cluster  area  to  obtain  the  imposed  vertical  velocity  wQ  in  (61) 
and  (62).  The  average  wQ  in  the  cloud  cluster  area  will  therefore  be 


The  atmospheric  sounding  taken  from  the  ship  Researcher  a 


1.7  cm  s  as  observed.  There  ls  no  horizontal  wind  initially.  Fig.  20, 
adopted  from  Ogura  et  al.  (1979),  shows  the  surface  wind  and  divergence  cal¬ 
culated  at  the  initiation  stage  of  the  cloud  band  at  00  GMT  12  August  197*1. 
It  is  evident  that  convective  cells  started  developing  in  the  area  of  the 
maximum  surface  convergence.  Almost  symmetric  decreases  of  surface  converg¬ 
ence  are  present  on  both  sides  of  the  cloud.  This  indicates  that  our  im¬ 
posed  circulation  wQ  described  above  is  fairly  justifiable. 

6.3.  A  24-hour  Integration 

In  this  section,  the  results  of  the  response  of  the  tropical  mesoscale 
cloud  band  to  the  large-scale,  low-level  lifting  process  through  the  appli¬ 
cation  of  the  present  cumulus  parameteri zation  scheme  will  be  illustrated  in 
terms  of  the  following  "circulation  descriptors" ,  as  defined  by  Kreitzberg 
and  Perky  (1977):  (1)  the  intensity  of  the  mesoscale  vertical  motion, 

(2)  the  life-cycle  of  the  development  of  the  mesoscale  circulation  and 

(3)  the  horizontal  width  scale  of  the  cloud  band  and  its  associated  precipi¬ 
tation  area.  For  verification  purposes,  the  model  results  will  be  compared 
with  the  observed  values  deduced  from  data  analysis. 

The  initial  conditions,  shown  in  Fig.  18,  are  conditionally  unstable 
except  for  a  stable  layer  around  700  mb.  There  is  no  initial  motion  field 
in  the  model  except  the  externally  imposed,  time- independent ,  large-scale, 
low-level  lifting  wQ  .  Since  the  model  is  symmetric  about  the  center  axis 
of  the  domain,  hereafter  all  of  the  model  results  are  presented  only  in  the 
right  half  of  tne  domain.  Also,  all  of  the  w  fields  presented  in  this 
section  do  not  include  the  imposed  large-scale  wQ  . 

During  the  first  six  hours  of  simulation,  it  is  observed  that  no  con¬ 
vective  activity  is  present,  primarily  because  the  large-scale  moisture 
convergence  has  not  reached  the  minimum  critical  value  for  convection  to 
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Fig.  20.  The  surface  wind  and  divergence  (in  the  units  of  10  :  s  ')  along 
with  radar  echoes  observed  at  00GMT  12  August  T 97u .  The  outer 
hexagon  Indicates  the  GATE  A/3  scale  ship  array  and  the  inner 
hexagon  the  3  scale  ship  array.  The  full  barb  for  wind  is 
5  m  s- ! .  (Adapted  from  Ogura  et  al.,  1979). 
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occur.  [Here,  the  minimum  critical  value  is  set  to  5  x  10  At  s  us  in 
Phillips  (1979).]  Fig.  21  shows  the  evolution  of  the  mesoscale  circulation 
in  t.ue  vector  field  of  [u,w]  (the  magnitude  of  w  has  been  enlarged  by  10 
times  because  of  the  scale  difference  between  u  and  w)  at  9,  13»  15,  and 
23  hrs .  At  9  hr,  the  flow  pattern  has  inflows  at  low  levels,  upward  motion 
at  the  central  portion  of  the  domain,  and  outflows  at  the  middle  level.  The 
upper  levels  have  almost  no  motion.  The  circulation  system  is  located 
within  160  km  of  the  center.  At  this  time,  we  can  categorize  this  mesoscale 
circulation  as  the  developing  stage.  As  the  convection  continues  to  de¬ 
velop,  the  magnitude  of  the  upward  motion  becomes  larger  and  the  upward 
motion  area  also  expands  and  extends  to  higher  levels.  At  13  and  15  hr, 
inflow  is  still  from  low  levels,  upward  motion  is  around  the  center  and  the 
outflow  is  at  upper  levels.  The  whole  circulation  system  is  much  larger 
than  at  the  developing  stage.  This  is  the  mature  stage.  Then  a3  the  con¬ 
vection  continues,  this  circulation  system  begins  to  split  and  the  intensity 
of  the  upward  motion  decreases  rapidly.  At  23  hr,  a  reverse  circulation  is 
present  at  the  low  levels,  categorizing  it  as  the  decaying  stage. 

The  reasons  for  the  above  evolution  of  the  mesoscale  circulation  were 
previously  explained  by  Yamasaki  (1984):  (1)  If  the  outflow  is  not  very 
strong  or  the  inflow  is  relatively  strong,  the  updraft  is  maintained  or 
intensified  at  almost  the  same  location;  (2)  if  the  outflow  of  the  downdraft 
air  is  strong  enough  for  a  new  updraft  to  occur  at  some  distance  from  the 
old  cloud,  the  cloud  has  a  short  lifetime  and  the  new  updraft  may  contribute 
to  tne  formation  of  another  cloud  system. 

Now  we  will  look  into  the  intensity  of  this  mesoscale  vertical  motion 
at  tne  developing,  mature,  and  decaying  stages.  Fig.  22  shows  the  evolution 
of  tne  vertical  velocity  at  different  stages.  It  is  seen  that  at  9  hr  the 
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Fig.  22.  Cros3  sections  of  the  predicted  vertical  velocity  at 
and  23  hours  with  a  contour  interval  of  2  cm/sec. 
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apward  motion  is  around  the  center.  Significant  downward  motion  with  a 
magnitude  of  about  2  to  3  cm/s  is  adjacent  to  an  upward  motion  area.  The 
level  of  maximum  vertical  velocity  is  located  at  around  2.7  km.  Above  the 
7.2  km  height,  the  upward  motion  is  close  to  zero.  At  this  stage  the  con¬ 
vection  is  dominated  by  middle  and  shallow  clouds.  This  is  due  to  the  foi — 
mulation  we  have  in  (43)  to  (45):  when  the  convection  is  weak  or  moderate, 
the  entrainment  effect  will  be  important;  therefore,  the  cloud  top  level  is 
not  located  at  a  higher  level.  At  13  and  15  hr,  the  intensity  of  upward 
motion  has  increased  and  both  the  upward  motion  area  and  the  downward  motion 
area  are  more  broad,  with  the  level  of  maximum  vertical  velocity  at  around 
10.3  km.  During  these  stages  the  convection  is  dominated  by  deep  clouds.  At 
23  hr,  two  features  are  noted:  one  is  the  downward  motion  appearing  in  the 
central  portion  of  the  domain,  while  another  is  the  secondary  maximum  of  the 
upward  motion  area  located  at  high  levels  some  distance  from  the  center.  As 
mentioned  earlier,  if  the  outflow  of  the  downdraft  air  is  strong  enough  for 
a  new  updraft  to  occur  at  some  distance  from  the  old  cloud,  the  new  updraft 
may  contribute  to  the  formation  of  another  cloud  system  with  this  secondary 
upward  motion  becoming  the  new  updraft.  This  secondary  upward  motion  con¬ 
sumes  most  of  the  moisture  supply  so  it  tends  to  reduce  the  original  convec¬ 
tion  around  the  center.  On  the  other  hand,  this  secondary  upward  motion 
area  cannot  develop  into  a  well-defined  convective  system  because  the  magni¬ 
tude  of  the  imposed  wQ  in  this  area  is  small.  Therefore,  the  low-level  mass 
convergence  is  insufficient  for  the  new  cloud  to  maintain  itself. 

Fig.  23  shows  the  predicted  cloud  heating  at  the  same  time  periods  as 
Fig.  22.  Basically,  the  cloud  heating  area  coincides  with  the  upward  motion 
area  in  Fig.  22.  At  9  hr,  the  heating  area  is  located  around  the  center  of 
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the  domain  and  stays  at  low  levels.  At  13  and  15  hrs,  the  level  of  maximum 
heating  has  shifted  to  high  levels  and  the  heating  area  has  become  more 
broad.  At  23  hr,  the  heating  area  has  moved  away  from  the  center. 

6.4.  Comparison  with  observation 

Fig.  24  shows  the  observed  vertical  velocity  (adopted  from  Ogura  et 
al.,  1979)  at  the  center  of  the  cloud  band  as  a  function  of  time  and  height, 
as  determined  by  wind  data  from  the  GATE  B-Scale  observational  network.  It 
is  evident  from  this  figure  that  prior  to  and  during  the  initiation  stage 
low-level  upward  motion  was  present  with  weak  subsidence  at  middle  and  upper 
levels.  As  convection  developed,  w  increased  in  intensity  and  the  area  of 
upward  motion  increased  greatly  in  vertical  extent.  The  local  maximum  w 
was  located  around  the  400-mb  level  during  the  mature  stage  of  the  cloud 
band  at  15  GMT.  The  decaying  stage  was  characterized  by  the  rapid  decrease 
in  upward  motion  at  low  levels  followed  by  the  development  of  descending 
motion,  while  significant  ascending  motion  remained  at  upper  levels. 

Fig.  25  shows  the  evolution  of  the  model-predicted  w  field  averaged 
horizontally  over  1 60  km  of  the  central  portion  of  the  domain,  which  is 
considered  the  width  of  the  cloud  band  simulated  by  the  model.  As  mentioned 
in  the  previous  section,  there  is  no  convection  present  initially.  The  con¬ 
vection  started  around  the  7th  hr  when  the  area-averaged  w  started  to 
increase  in  intensity.  The  level  of  maximum  w  remained  at  the  2.7  km 
level  until  12  hr.  After  12  hr,  the  area-averaged  w  rapidly  increased  its 
magnitude  and  the  level  of  maximum  w  shifted  to  a  higher  altitude.  At  16 
hr,  the  upward  velocity  reached  its  maximum  value  of  around  10  cm/sec  at  9.9 
km  height.  Around  i3  hr,  the  averaged  w  rapidly  decreased  in  intensity  and 
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ig.  25.  Tiae-heignt  sec-ion  cf  the  model-predicted  vertical  velocity  aver¬ 
aged  over  the  central  1 50  km  of  the  model  domain  with  a  contour 
interval  of  0.5  cm/sec. 


low-level  downward  motion  was  observed.  Then  the  whole  system  gradually 


started  decaying.  At  24  hr,  the  lower  half  of  the  domain  was  dominated  by 
downward  motion. 

Comparing  the  simulated  w  field  with  observations,  one  can  see  that 
the  model  results  agree  well  with  the  observations  in  that  (1)  the  area- 
averaged  w  increases  in  intensity  with  time,  (2)  prior  to  the  mature  stage 
the  level  of  maximum  w  is  located  at  low  levels;  then  as  the  convection 
developed  the  location  of  maximum  w  is  shifted  to  a  higher  altitude,  i.e., 
the  height  of  the  local  maximum  of  w  increases  with  time,  and  (3)  after 
the  mature  stage  the  intensify  of  upward  motion  at  low  levels  decreases 
rapidly  and  then  becomes  downward  motion. 

Fig.  26  shows  the  model-predicted  surface  precipitation  as  a  function 
of  time  and  distance  from  the  center.  The  width  and  distribution  of  pre¬ 
cipitation  almost  coincide  with  those  of  the  predicted  upward  motion.  Fig. 
27  shows  the  observed  precipitation  accumulation  for  the  cloud  band  from  12 
to  18  GMT  12  August  1974.  The  observed  average  rainfall  accumulation  in  the 
center  of  the  cloud  band  (8.5  to  10  N  and  23.5  to  22  W)  is  around  11.8  mm 
from  12  to  13  GMT.  The  area-averaged  rainfall  from  our  simulation  during 
the  same  period  is  11.9  mm,  which  is  very  close  to  the  observation. 

For  comparison,  Fig.  28  shows  the  time-height  section  of  the  the  pre¬ 
dicted  vertical  velocity  averaged  over  the  central  160  km  of  the  model  do¬ 
main  with  the  use  of  a  combination  of  Kuo's  1  965  and  1  974  schemes  (see 
Section  3)  with  the  value  of  b  set  to  zero.  A  comparison  of  this  figure 
with  Fig.  20  shows  that  as  the  convection  started  the  level  of  maximum  up¬ 
ward  motion  immediately  shifted  to  i  higher  level  and  then  remained  at  the 
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Model-predicted  rainfall  accumulation  as  a  function  of  time  and 
distance  from  the  center  in  units  of  mm/hr  with  a  contour  interval 
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The  observed  rainfall  accumulation  for  the  period  from  1 2GM' 
7  8G'-^T  12  August  1974  in  units  of  mm. 
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the  intensity  of  w  decreased  rapidly;  however,  no  downward  motion  was  ob¬ 
served.  At  23  hr,  the  averaged  w  started  to  increase  its  intensity  in¬ 
stead  of  entering  the  decaying  stage.  In  this  simulation,  one  can  see  that 
(1)  the  height  of  the  local  maximum  w  did  not  increase  with  time,  and  (2) 
the  downward  motion  did  not  appear  at  low  levels  during  the  decaying  stage. 

7.  Concluding  Remarks 

The  determination  of  the  vertical  distributions  of  heating  and  mois¬ 
tening  bv  the  cumulus-scale  motions  and  the  estimation  of  reasonable  rain¬ 
fall  rates  are  some  of  the  major  problems  related  to  the  parameter i zat ion  of 
cumulus  convection.  Kuo  has  proposed  that  both  the  amount  and  the  vertical 
distribution  of  convective  heating  by  convergence-controlled  cumulus  convec¬ 
tion  can  be  calculated  by  using  the  moisture  convergence  Mt  and  the  tempera¬ 
ture  difference  (T  -  T)  between  the  cloud  and  the  environment.  In  this 
report,  a  combination  of  Kuo's  1965  and  1 97 ^  scheme?  were  first  applied 
using  a  semiprognostic  approach  without  considering  the  entrainment  to  the 
GATE  Phase  III  data.  It  was  found  that  the  vertical  heating  profiles  are  in 
good  agreement  with  observations  during  periods  of  strong  convective  condi¬ 
tions.  However,  they  are  in  poor  agreement  during  periods  of  weak  convec¬ 
tion.  This  problem  can  be  resolved  by  considering  the  vertical  extent  of 
the  cumulus  cloud  (entrainment  effect)  and  the  vertical  profiles  of  heating 
based  on  different  convective  conditions.  This  led  to  the  present  scheme. 
Comparisons  between  the  present  scheme  and  Kuo's  scheme  without  entrainment 
were  then  made  using  semi  prognostic  tests,  Mowing  that  the  present  scheme 
prose 0*3  a  drastic  improvement  in  cloud  heating  profiles  during  periods  of 
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weak  convection.  Generally,  the  present  scheme  can  reproduce  the  observed 
rainfall  and  heating-drying  profiles  quite  well  in  the  semi  prognostic  tests 
using  the  GATE  Phase  III  data. 

A  prognostic  application  of  the  present  modified  Kuo  scheme  was  made  by 
incorporating  the  scheme  into  a  two-dimensional,  cloud  cluster  model  to 
simulate  a  tropical  cloud  band  under  the  influence  of  a  low-level,  large- 
scale  lifting.  A  realistic  life  cycle  behavior  was  simulated  and  the  model 
results  from  our  present  scheme  bear  considerable  resemblance  to  many 
aspects  of  the  observations  including  (1)  the  evolution  of  an  area-averaged 
vertical  velocity,  (2)  the  intensity  and  distribution  of  precipitation,  and 
(3)  the  development  of  low-level  downward  motion  during  the  decaying  state. 
The  results  indicate  that  (1)  linking  rainfall  to  a  fraction  of  the  large- 
scale  moisture  supply  is  reasonable  in  mesoscale  tropical  cloud  band  simula¬ 
tion,  and  (2)  the  incorporation  of  the  entrainment  effect  into  Kuo’s 
approach  appears  to  be  essential  to  reproduce  the  evolution  of  the  area- 
averaged  vertical  velocity  and  the  life  cycle  of  the  convective  system. 
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PART  II 

An  Evaluation  of  the  Cumulus  Parameterization  Scheme  using  the  Air  Force 
Geophysics  Laboratory  Global  Moist  Forecast  Model 

1 .  Introduction 

A  good  cumulus  parameterization  scheme  should  be  based  on  understanding 
cloud-environment  interactions.  These  interactions  include  (a)  the  control¬ 
ling  influence  of  environmental  conditions  on  the  development  and  intensity 
of  an  ensemble  of  cumulus  clouds  and  (b)  the  physical  processes  through 
which  this  cumulus  ensemble  modifies  its  environment.  The  verification  of  a 
cumulus  parameterization  scheme  should  also  be  based  on  its  ability  to  real¬ 
istically  represent  the  heating  and  drying  effects  of  cumulus  clouds  and 
predict  accurately  the  cumulus  precipitation.  In  Scientific  Report  Mo.  1 
(Kao  and  Ogura,  1985),  we  evaluated  the  Arakawa-Schubert  (A-S)  scheme  of 
cumulus  parameterization.  In  Part  One  of  this  report,  we  evaluated  the  Kuo 
scheme.  Two  methods  were  used  in  these  evaluation  processes.  One  is  the 
semi-prognostic  approach,  in  which  the  parameterized  cumulus  heating  and 
drying  profiles  and  precipitation  are  compared  with  the  values 
diagnostically  determined  from  the  heat  and  moisture  budget  equations.  The 
second  approach  is  to  incorporate  the  cumulus  parameter i zat ion  schemes  into 
a  fully  prognostic  cloud  cluster  model.  The  ability  of  cumulus  parameteri¬ 
zation  schemes  to  simulate  the  cloud-environment  interaction  processes  and 
to  reproduce  the  life  cycle  of  a  cloud  cluster  can  be  evaluated.  Through 
these  studies,  an  improved  algorithm  to  implement  the  A-S  scheme  and  an 
improved  Kuo's  scheme  including  the  entrainment  effect  were  developed. 

The  ultimate  test  of  any  cumulus  par ameter i zat i on  scheme  is  to  actually 
incorporate  that  scheme  into  a  large-scale  prediction  model.  The  advantage 
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of  this  approach  is  obviously  that  the  predicted  cloud  coverage  and  precipi¬ 
tation  amount  can  be  compared  with  actual  observations.  The  drawback  is 
that  when  the  prediction  is  incorrect,  it  is  not  certain  whether  the  fault 
rests  on  a  poor  cumulus  parameterization  scheme  or  on  the  large-scale  pre¬ 
diction  model  itself.  In  spite  of  these  shortcomings,  this  test  is  neces¬ 
sary  for  the  improvement  of  numerical  weather  prediction.  This  part  of  the 
report  presents  the  results  of  applying  the  Kuo  scheme  and  the  A-S  scheme  of 
cumulus  parameterization  to  the  AFGL  Global  Moist  Forecast  Model.  Section  2 
provides  a  brief  description  of  the  AFGL  model.  Section  3  presents  the 
prediction  results,  followed  by  the  discussion  of  the  results  in  Section  4. 
2.  The  AFGL  Global  Moist  Forecast  Model 

The  equations  and  framework  of  the  AFGL  Global  Moist  Forecast  Model  are 
described  in  Brenner  et  al .  (1982).  A  brief  description  will  be  provided 

here  as  background.  The  AFGL  Global  Moi3t  Forecast  Model  is  a  spectral 
model  with  12  layers  and  rhomboidal  30  waves.  The  prognostic  equations  are 
vorticity,  divergence,  temperature,  specific  humidity,  and  the  logarithm  of 
the  surface  pressure  (continuity  equation).  The  vorticity  and  divergence 
are  given  conventionally  as  Laplacians  of  a  streamf unct ion  and  velocity 
potential  respectively.  The  velocity  components  and  geopotential  are  ob¬ 
tained  diagnostically  with  the  aid  of  the  hydrostatic  relation.  The  verti¬ 
cal  coordinate  of  the  model  is  o  =  1  -P/p*  .  defined  by  Phillips  (1959), 
with  a  specification  of  layer  locations  described  by  Brown  (1974)  and  Phil¬ 
lips  (1975).  The  numerical  methods  in  the  model  include  spectral  represen¬ 
tation  in  the  horizontal,  the  Arakawa  quadratic  conserving  finite 
differencing  in  the  vertical,  and  the  semi- impl ici t  time  integration  scheme. 
The  physical  effects  included  in  the  model  are  the  influences  of  orography, 
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position-dependent  surface  friction,  moisture  physics,  and  subscale  horizon¬ 
tal  dissipation,  parameterized  by  diffusion.  Evaporation  and  sensible  heat 
flux  from  the  oceans  are  also  included.  The  application  of  the  moisture 
physics  consists  of  a  sequence  of  three  steps  to  adjust  the  temperature  and 
specific  humidity.  Each  step  possesses  a  characteristic  spatial  scale:  (1) 
cumulus  convection  in  conditionally  unstable,  generally  unsaturated  large- 
scale  flow,  (2)  large-scale  condensation  in  stable,  saturated  large-scale 
flow,  and  (3)  dry  convection  in  unstable,  unsaturated  large-scale  flow.  The 
adjustments  are  applied  after  the  provisional  values  of  the  dynamic  varia¬ 
bles  at  a  new  time  step  are  obtained  through  time  integration.  Large-scale 
condensation  produces  nonconvecti ve  and  stratiform  precipitation.  The  ef¬ 
fects  of  the  cumulus  clouds  are  parameteri zed  using  Kuo's  modified  1965 
scheme.  The  highlights  of  this  scheme  consist  of  computing  vertically  inte¬ 
grated  moisture  convergence  in  the  seven  lowest  layers  of  each  column.  If 
the  computed  amounts  exceed  10  '  s  ,  and  the  atmosphere  is  unstable,  then 
cumulus  clouds  will  be  allowed  to  develop.  An  unstable  region  is  defined  as 
extending  from  the  bottom  to  the  first  layer,  which  is  warmer  than  a  moist 
adiabatically  lifted  parcel.  If  such  an  unstable  column  is  found,  the  fol¬ 
lowing  sums  are  computed  from  the  bottom  layer  to  the  top  of  the  unstable 
part  of  the  column: 
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where  the  term  "WATER"  in  the  expression  for  Q  ^  is  the  moisture  converg¬ 
ence  in  the  unstable  column,  qc  and  Tc  are  the  cloud  mixing  ration  and  tem¬ 
perature  calculated  from  the  lifted  parcel  values,  and  L  is  the  latent  heat. 
DTKUOk  represents  the  latent  heat  release  of  convective  precipitation  and  is 
added  to  the  temperature  forecast  before  cumulus  parameter! zation.  DQKUOk 
represents  the  change  in  specific  humidity  by  cumulus  parameterization  and 
is  added  to  the  value  at  the  previous  time  step.  In  this  version  of  the  Kuo 
scheme,  the  condensed  water  is  also  allowed  to  evaporate  while  falling  into 
the  lower  layers.  Fig.  29  contains  the  flow  chart  of  Kuo's  cumulus  parame¬ 
terization  scheme  in  the  AFGL  spectral  model. 

3.  Results 

The  experiments  consist  of  3  cases,  each  case  consisting  of  4  integra¬ 
tions  of  the  AFGL  spectral  model  for  24-hour  simulations  and  starting  from 
identical  conditions.  Integration  1  uses  Kuo's  parameterization  in  the 
original  AFGL  spectral  model,  as  described  above.  Integration  2  uses  the 
same  Kuo  scheme,  but  the  condensed  water  i3  not  allowed  to  evaporate.  Inte¬ 
gration  3  uses  the  present  modified  Kuo  scheme  reported  in  Part  One  of  the 
report,  and  Integration  4  uses  the  A-S  ( Arakawa- Schubert )  cumulus  parameter¬ 
ization  scheme.  A  detailed  description  of  the  A-S  scheme  is  given  in  Scien¬ 
tific  Report  No.  1.  The  initial  data  used  for  the  global  test  forecasts 
consisted  of  FGGE  Level  1 1 1— 3  on  1,  4,  and  9  April  1979.  Since  the  rainfall 
distribution  is  the  major  feature  that  needs  to  be  simulated  accurately  by  a 
"good"  formulation  of  parameterization,  only  the  rainfall  prediction  will  be 
addressed  in  this  report.  In  view  of  the  limited  observed  precipitation 
data  available,  the  discussion  of  forecast  results  will  be  limited  to  the 
area  of  the  U.S.  where  the  observed  hourly  precipitation  data  are  available 
through  the  National  Center  for  Atmospheric  Research. 
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KUO  (1965)-Type  Convection 


condition: 


Fig.  29.  The  schematic  diagram  of  the  computational  procedure  of  the  Kuo-type 


(  1965)  scherm.  in  the  AFC!-  model. 


3.1.  The  April  1,  1979  case 


Starting  with  an  initial  condition  from  the  FGGE  Level  III-B  data  at  00 
GMT  April  1,  1979,  a  2^-hour  simulation  was  undertaken  with  four  different 
versions  of  cumulus  parameterization.  In  the  first  few  hours  of  the  simula¬ 
tion,  any  adjustments  to  the  initial  data  may  cause  unrealistic  rainfall 
patterns.  Therefore,  the  comparison  between  the  predicted  rainfall  and 
observations  will  be  made  only  for  the  second  12  hours. 

Fig.  30  shows  the  surface  synoptic  conditions  at  1200  GMT  1  April  1979. 
The  major  features  included  a  low-pressure  system  centered  around  the  Texas- 
Oklahoma  border.  A  cold  front  associated  with  this  low  was  located  at  the 
southern  end  of  this  low  center  and  extended  to  Texas.  A  stationary  front 
was  located  at  the  northeastern  end  of  this  low  center  and  extended  to  the 
East  Coast.  Also  there  was  a  cold  front  lying  across  North  Dakota,  South 
Dakota,  Nebraska,  and  Wyoming.  During  the  time  between  1200  and  2^00  GMT., 
the  low  pressure  system  moved  northeastward  to  the  Arkansas-Missouri  border, 
while  the  associated  cold  front  followed  the  movement  of  this  low  and  gradu¬ 
ally  became  a  stationary  front.  Also,  a  squall  line  system  developed  to  the 
east.  Heavy  convective  precipitation  occurred  with  this  frontal  system  from 
Kentucky  southwest  to  Arkansas  and  Louisiana  during  this  12-hr  period,  as 
shown  in  the  radar  summary  in  Fig.  31  •  The  stationary  front  northeast  of 
this  low  center  also  moved  along  with  this  low  and  produced  moderate  rain¬ 
fall  from  Massachusetts  southwest  to  Kentucky.  The  cold  front  around  Ne¬ 
braska  and  Wyoming  contributed  some  light  rainfall  over  Montana,  Wyoming, 
Colorado,  Kansas,  and  Nebraska. 

Fig.  32  shows  the  observed  12-hr  accumulated  rainfall  ending  at  2^  GMT 
1  April  1979  over  the  U.S.  NO  data  is  available  over  the  oceans  and  Canada; 
therefore,  the  contour  lines  over  these  areas  are  produced  by  interpolation 
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and  are  not  meaningful..  This  figure  shows  a  large  amount  of  accumulated 
rainfall  from  Massachusetts  southwestward  to  Texas  and  some  light  rainfall 
over  Colorado,  Wyoming,  and  Montana.  The  rainfall  areas  coincide  with  the 
areas  of  radar  echoes  in  Fig.  31 . 

Fig.  33  is  the  predicted  12-hr  accumulated  total  rainfall  from  Integra¬ 
tion  1.  The  display  area  is  chosen  to  be  the  same  as  in  Fig.  32.  The  fore¬ 
cast  total  precipitation  includes  the  convective  and  large-scale 
precipitation.  In  this  integration,  the  forecast  rainfall  area  is  from  Mas¬ 
sachusetts  west  to  Wyoming  and  Colorado  with  heavy  precipitation  over  Indi¬ 
ana  and  Illinois.  A  small  area  of  light  rainfall  is  predicted  over 
Washington  and  Oregon.  This  predicted  light  rainfall  area  agrees  in  general 
with  the  observation;  but  heavy  precipitation  is  not  predicted  over  Arkansas 
and  Louisiana,  as  observed.  The  predicted  rainfall  in  this  integration  is 
almost  totally  contributed  by  large-scale  condensation.  The  convective 
rainfall  (not  shown)  is  very  small,  i.e.,  not  over  0.1  mm  during  this  pe¬ 
riod.  Fig.  3^  is  the  predicted  rainfall  from  Integration  2.  Fig.  3^a  is  the 
predicted  total  precipitation.  The  forecast  rainfall  area  is  essentially 
the  same  as  thac  from  Integration  1  except  the  rainfall  area  is  slightly 
bigger.  The  predicted  convective  rainfall  (Fig.  3^b),  is  located  from  the 
East  Coast  to  Missouri  and  Arkansas  in  the  broad  area  of  observed  heavy 
precipitation.  However,  the  amount  is  still  very  small  compared  with  that 
from  the  large  scale  condensation.  Fig.  35  is  from  Integration  3.  Compar¬ 
ing  Fig.  35a  with  Fig.  32,  one  can  see  that  the  predicted  rainfall  area,  es¬ 
pecially  the  rainfall  associated  with  the  cold  front,  agrees  much  better 
with  the  observation.  Fig.  35b  shows  that  the  rainfall  area  from  Kentucky 
to  Louisiana  is  produced  purely  by  convective  rainfall.  This  convective 
rainfall  area  almost  coincides  with  the  intense  radar  echo  area  in  Fig.  31; 
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Fig.  3a.  Model-predicted  12- hr  rainfall  ending  /NOO  GMT  1  April  1979 
interpolated  over  the  U.S.  from  Integration  2  with  a  contour  in¬ 
terval  of  5  mm:  (a'1  total  precipitation,  (b)  convective  rainfall. 
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however,  the  amount  of  the  predicted  rainfall  rate  is  still  smaller  than  the 
observation  and  the  center  of  maximum  rainfall  is  too  far  to  the  northeast. 
Fig.  36  shows  the  results  from  Integration  4.  In  Fig.  36a,  the  predicted 
total  rainfall,  the  rainfall  area  covers  most  of  the  U.S.  except  for  the 
West.  The  maximum  rainfall  area  is  over  the  Minnesota-Iowa  border,  which  is 
too  far  northwest  compared  with  the  observation.  This  rainfall  is  mostly 
contributed  by  cumulus  clouds,  as  parameterized  by  the  A~S  scheme  (Fig. 
36b).  Another  area  of  convective  precipitation  is  predicted  in  a  band  from 
Tennessee  to  Louisiana  close  to  the  observed  rainfall  area. 

3.2.  The  April  4,  1979  case 

The  second  24-hr  forecast  discussed  in  this  report  begins  at  00  GMT  4 
April  1979.  Fig.  37  shows  the  surface  synoptic  chart  at  1200  GMT  4  April 
1979.  The  main  features  at  this  time  included  a  low-pressure  system  cen¬ 
tered  over  northwest  Kentucky,  a  stationary  front  over  Alabama  and  a  cold 
front  over  the  North  Dakota-Canadian  border  which  extended  west  and  north¬ 
westward  to  North  Dakota,  Montana,  and  into  Canada.  During  the  next  12  hr, 
the  low-pressure  system  moved  northeastward  to  the  Ohio-Pennsyl vania  border 
and  produced  non-convecti ve  rainfall  over  the  East  Coast.  The  stationary 
front  over  Alabama  also  slowly  moved  northeastward  to  Georgia  with  the 
southern  part  of  this  front  gradually  becoming  a  cold  front  during  this 
period.  An  intensive  radar  echo  associated  with  this  cold  front  extended 
the  full  length  of  the  Atlantic  Coastal  States  (Fig.  38).  The  cold  frontal 
system  over  the  North  Dakota-Canadian  border  had  also  moved  nor theastward 
into  Canada.  Several  small  radar  echoes  associated  with  this  front  were 
located  over  North  Dakota,  South  Dakota,  Wyoming,  and  Montana. 
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Fig.  39  shows  the  observed  12-hr  accumulated  rainfall  ending  at  2400 
GMT  4  April  1979  over  the  U.S.  Heavy  rainfall  over  the  East  Coast  was  again 
evident  with  the  heaviest  area  of  precipitation  over  northern  Florida. 
Light  rainfall  existed  from  the  Iowa-Missouri  border  northwest  to  Washington 
State.  Also  light  rainfall  was  present  over  Texas. 

Fig.  40  is  the  predicted  12-hr  accumulated  total  rainfall  over  the  U.S. 
from  Integration  1.  The  forecast  precipitation  area  is  concentrated  over 
the  northeast  with  light  rainfall  being  predicted  over  the  northwestern  U.S. 
The  precipitation  areas  agree  somewhat  with  the  observations  (Fig.  39), 
except  it  totally  missed  the  heavy  rainfall  over  Florida  and  Texas.  The 
predicted  rainfall  is  mostly  due  to  large-scale  condensation  and  the  convec¬ 
tive  rainfall  from  this  integration  is  again  very  small.  Fig.  41  is  the 
predicted  rainfall  from  Integration  2  and  Fig.  41a  is  the  total  rainfall. 
The  predicted  rainfall  area  is  essentially  the  3ame  as  that  from  Integration 
1  except  it  expanded  in  area  from  the  northeast  down  to  Florida.  In  Fig. 
41b,  predicted  convective  rainfall  over  Florida  is  now  evident.  Also, 
light  convective  rainfall  is  predicted  over  the  borders  of  Montana,  North 
Dakota,  and  Canada.  However,  the  amount  of  convective  rainfall  is  very 
small  compared  with  that  from  lar  , e-scale  condensation.  For  Integration  3 
(Fig.  42),  the  predicted  total  rainfall  has  two  centers  over  the  East  Coast 
and  is  in  good  agreement  with  thr  observation  (Fig.  "'9).  Fig.  42b  clearly 
shows  that  the  areas  of  rainfall  over  riorida  and  Texas  are  the  result  of 
convection  only.  In  Fig.  43,  from  Integration  4,  the  precipitation  area  is 
concentrated  over  the  Northeast  Coast  and  over  the  ocean  off  of  Florida. 
Some  light  rainfall  is  predicted  over  Washington  State  and  extending  ea3t 
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Model-predicted  12-hr  rainfall  ending  at  2403  GMT  4  April  1 979 
interpolated  over  the  U.S.  from  Integration  2  with  a  contour  in¬ 
terval  of  5  mm:  (a)  total  rainfall,  (b)  convective  rainfall. 
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Fig.  42 .  As  in  Fig.  41  except  for  Integration  3- 
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and  southeast  to  Texas.  It  also  shows  some  rainfall  over  the  West  Coast. 
The  forecast  rainfall  area  in  general  agrees  with  the  observations  but  cov¬ 
ers  an  area  much  larger  than  the  observation. 

3.3.  The  April  9,  1 979  case 

The  third  24-hr  forecast  presented  in  this  report  starts  at  00  GMT  9 
April  1979.  Fig.  44  shows  the  surface  synoptic  conditions  at  1200  GMT  9 
April  1979.  The  major  features  included  a  low-pressure  system  centered  over 
West  Virginia  and  a  cold  front  associated  with  this  low  located  to  the  south 
of  this  low  center  and  extending  southward.  Another  cold  front  was  located 
over  the  southwestern  part  of  Canada  and  extended  south  and  southwestward  to 
California.  During  the  time  between  1200  and  2400  GMT,  the  low-pressure 
system  moved  northeastward  and  the  radar  echo  associated  with  this  system 
was  over  the  East  Coast  (Fig.  45).  The  cold  front  over  Canada  moved  south- 
|  ward  to  Montana  and  southwestward  to  Wyoming,  Utah,  Nevada,  and  California. 

Fig.  46  shows  the  observed  12-hr  accumulated  rainfall  ending  at  24  GMT 
9  April  1979.  Precipitation  covers  an  area  from  Washington  State  southeast- 
|  ward  down  into  Texas  and  Louisiana.  Heavy  rainfall  was  produced  over  the 

|  N:r" re ast . 

Fig.  *«'’  is  the  predicted  12-hr  accumulated  total  rainfall  from  Integra- 
)  1 1  or.  '.  The  predicted  rainfall  areas  are  mainly  over  the  northeast  with 

‘  p  it  r.y  ire  is  ?f  1  i grit  rainfall  extending  from  Washington  State  southeastward 

.  to  T-'X.is.  As  before,  the  convective  rainfall  from  this  integration  is  very 

► 

P  small.  Dveral 1 ,  the  p^eci pi  tat  ion  pa* tern  matches  that  in  Fig.  46  reasona- 

>  Ply  well.  Fig.  43  shows  the  forecast  rainfall  from  Integration  2.  Due  to 

> 

j.  con  veer  i  on,  tr.e  predicted  rainfall  areas  over  the  northeast  expanded  south- 

I  ward  tve-  the  Atlantic  coastal  waters.  The  rainfall  area  over  the  central 


United  Ctates  was  also  mainly  due  to  convection. 


Fig.  49,  from  Integration 
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Fig.  47. 


Model-predicted  12-hr  rainfall  ending  at  2400  GMT  9  April  1979 
interpolated  over  the  U.S.  from  Integration  1  with  a  contour  in¬ 
terval  of  5  mm. 
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Fig.  ^9.  As  in  Fig.  48  except  for  Integration  3. 
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3,  shows  that  heavy  rainfall  was  produced  over  the  Atlantic  coastal  area. 
Any  comparison  with  observations  over  that  area  would  be  misleading  because 
no  observations  were  available.  Fig.  50  shows  the  forecast  rainfall  from 
Integration  *1.  In  general,  precipitation  covers  the  East  Coast  and  the  High 
Plains  area.  The  area  of  maximum  rainfall,  over  the  Wyoming-Colorado  border, 
has  a  value  above  30  mm,  which  is  much  too  high  compared  with  the  observa¬ 
tion.  This  rainfall  is  contributed  totally  by  convection  due  to  the  A-S 
scheme . 

4.  Discussion 

The  results  in  the  preceding  sections  have  led  us  to  the  following 
conclusions  from  the  four  different  cumulus  parameterization  schemes: 

(a)  In  the  original  AFGL  spectral  model,  Kuo's  1965  scheme  was  used. 
In  addition,  the  condensed  falling  water  was  allowed  to  evaporate  into  the 
lower  layers  if  they  were  not  saturated.  This  version  of  the  scheme  pro¬ 
duced  almost  no  convective  rainfall  and  the  total  rainfall  was  mainly  con¬ 
tributed  by  the  large-scale  condensation.  The  lack  of  convective  rainfall 
was  caused  by  the  rainfall  evaporating  before  it  reached  the  ground. 

(b)  Kuo’s  1965  scheme  without  any  evaporation  of  falling  rain  produced 
some  convective  rainfall.  However,  the  amount  was  still  very  small  compared 
to  that  from  the  large-scale  condensation.  As  discussed  in  Part  One  of  the 
report,  this  version  of  the  scheme  used  a  large  portion  of  the  large-scale 
supply  of  moisture  to  moisten  the  atmosphere,  leaving  a  small  portion  to 
precipitat ion. 

(c)  The  present  modified  Kuo  scheme  produced  a  reasonable  amount  of 
convective  rainfall  because  the  'b'  parameter  in  this  scheme  is  always  very 
small.  The  predicted  rainfall  area  agreed  very  well  with  the  observations 
in  all  three  of  the  test  cases. 
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(d)  The  A-S  scheme  tended  to  produce  more  convective  rainfall  than 
other  cumulus  parameterization  schemes.  The  predicted  rainfall  area  from 
this  scheme  was  larger  than  the  observations  in  the  above  3  cases. 

This  work  clearly  demonstrated  the  ability  of  the  AFGL  model  to  predict 
large-scale  stratiform  precipitation.  With  the  incorporation  of  the  present 
modified  Kuo  scheme,  the  area  of  convective  precipitation  can  also  be  well 
predicted.  However,  the  predicted  convective  precipitation  area  is  genei — 
ally  broader  and  the  amount  of  rainfall  smaller  than  observations.  This 
problem  may  be  attributed  to  the  resolution  of  the  AFGL  model  in  that  it 
cannot  resolve  the  narrow  band  of  the  cold  front,  the  major  mesoscale  rain- 
producing  system.  The  implementation  of  the  A-S  scheme  in  the  AFGL  model 
also  produced  a  reasonable  distribution  of  convective  precipitation.  The 
precipitation  area  is  more  concentrated  and  is  sometimes  produced  in  an 
observed  clear  area.  This  result  may  have  been  caused  by  the  current  imple¬ 
mentation  of  the  A-S  scheme  in  the  AFGL  model,  in  which  the  cloud  base  is 
assumed  to  be  500  m  above  the  ground  surface.  This  assumption  is  valid  in 
the  tropics,  but  it  may  produce  unrealistic  convective  rainfall  in  dry 
areas . 

Future  improvements  in  precipitation  forecasting  using  the  AFGL  model 
include  the  prediction  of  mesoscale  dynamics  related  to  strong  convection 
and  precipitation.  A  major  task  in  this  direction  is  to  incorporate  the 
prediction  of  the  diurnal  temperature  variation  which  can  significantly 
affect  the  structure  of  mesoscale  systems  and  the  time  and  location  of  pre¬ 
cipitation.  The  inclusion  of  boundary  layer  prediction  and  radiation,  which 
have  already  been  undertaken  by  AFGL  concurrently  with  this  project,  should 
fulfill  this  requirement.  Once  all  of  these  physical  processes  are  inte¬ 
grated  to  the  AFGL  model,  another  task  to  improve  the  prediction  is  to 
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increase  the  model  resolution.  This  task  can  be  achieved  by  either  increas¬ 
ing  the  wave  number  retained  in  the  AFGL  model  or  using  a  nested  technique 
for  limited  area,  fine  grid  prediction.  A  more  comprehensive  test  of  the 
fine  grid  model  with  an  extended  period  of  prediction  will  therefore  be 
needed. 
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